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Abstract 

In this paper, we introduce a modified epidemic model on regular and scale-free 
networks respectively. We consider the birth rate 5, cure rate 7, infection rate A, 
a from the infectious disease, and death rate (3 from other factors. Through mean- 
field analysis, we find that on regular network there is an epidemic threshold A c 
dependent on the parameters 5, 7, a, and j3; while for power law degree distribution 
network epidemic threshold is absent in the thermodynamic limit. The result is 
the same as that of the standard SIS model. This reminds us the structure of the 
networks plays a very important role in the spreading property of the infectious 
disease. 
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1 Introduction 

During the past few years, the spread of disease has been one of the focuses in 
the field of statistical physics. A a great deal of epidemiological research work 
has been done on various networks. Two epidemic models SIS and SIR have 
been widely studied[l-9]. In these models, each node of networks represents 
an individual and each link is the connection along which the individuals 
interact and the disease can spread. For SIS epidemic model, each individual 
can exist in two possible states: susceptible (or healthy) and infected. At each 
time step, each healthy individual can be infected at rate A if there is one 
or more infected individuals in its nearest neighbors. Meanwhile, an infected 
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individual may recover and become susceptible at rate 7. The SIR model 
assumes that individuals can exist in three possible states: susceptible (or 
healthy), infected and removed. The main difference from the SIS model is 
that once an individual gets infected, it is removed and can not be infected 
any more. It is easy to understand that both the properties of disease and 
network topology determine the dynamical behaviors of the disease spreading. 
Studies of SIS model and SIR model show that, on regular networks there is 
an epidemic threshold A c . If the effective spreading rate A > A c , the infection 
spreads and becomes endemic. Otherwise, the infection will disappear. While 
the epidemic threshold is absent on scale-free networks in the thermodynamic 
limit[10-12]. 



From the definitions of SIS and SIR models, we know that both SIS and SIR 
models assume the number of individuals to be a constant. However, some 
disease may cause enough deaths to influence the population size. So it is 
necessary to take the birth and death rates into account. In this paper, we 
introduce a modified model, considering the birth rate 5, treatment rate 7, 
infection rate A, and two death rates: j3 due to this infectious disease and a 
due to other factors, which will be described later. Our work is to study the 
influence of above parameters to the epidemic thresholds on different complex 
networks. 



2 Model 



We think of our individuals as being spatially distributed on the network Z. 
Each site of Z is empty or occupied by at most one individual. We give each 
site a number: 0, 1 or 2. They describe empty state, a healthy individual 
occupation and an infected individual occupation respectively. The state of 
the system at time t can be described by a set of numbers, 0, 1, 2. That means 
if the system is in state A and the site x G Z, then A t (x) G {0, 1,2}. Each 
site can change its state with a certain rate. An empty site can give birth to a 
healthy individual at rate 5. A healthy individual can be infected by contact 
at rate A if there are infected individuals in its nearest neighbors, or die at rate 
a due to other factors. An infected individual can be cured at rate 7 or die at 
rate (5 due to this infectious disease. If an individual dies, there is an empty 
site left. Of course, each site can also maintain its state. We define rii(x,t) as 



2 



the number of the nearest neighbors of site x in state i at time t. 



— > 1 at rate 5 

1 — > at rate a 

1 — > 2 at rate n 2 A 

2 — > 1 at rate 7 
2 -> at rate /3 



In above expressions, 5, a, [3, 7 and A are all non-negative. We assume a is 
relatively very small. The expression n 2 A means that a healthy individual with 
n 2 infected nearest neighbors gets infected at rate n 2 A. Not difficult to see, 
if 5, a and (3 equal 0, this model turns to SIS model; if 5, a, and 7 equal 0, 
this model turns to SIR model. If a and 7 equal 0, this model turns to "forest 
fire", which has been widely studied[13]. 



2.1 Mean-field Method on regular networks 



First, we solve the model by mean-field method on regular network without the 
consideration of spatial fluctuation. We use the density x and y (x,y G [0, 1]) 
to replace the numbers of the healthy individuals and the infected individuals 
respectively. "n 2 A" can be replaced as "X(k)y , \ where (k) is the average 
number of the nearest neighbors of one node. The evolution equations of x 
and y are governed by: 

Ox 

— = [l - x - y )S - ax - X(k)xy + ^y (1) 
dy 

— = A (k) xy - 7y - f3y (2) 



In Eq.(l), the expression (1 — x — y) is the density of empty site, (k) y is the 
probability that the nearest neighbors of one healthy individual are infectious. 

Let || = and |jjf = 0, we get the steady-state solutions: 
(I) 

5 

x = — — =, y = 0; (3) 



a + 5' 



and 
(II) 



2±£ £A(fc)-(£ + a)(7 + /3) 
X A (k) ' V \{k)(6 + P) [ ' 
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Now, I will do stability analysis. For solution (I), the Jacobean is: 

(5) 



J = 



(-a-5 i-5- s -^ 

i 5+a 



\ ^-(7 + /3) 
The determinant and the trace of J: 



|J| = -(« + *)[^-(7 + fl] (6) 
o + a 



Tr(J) = -a-5+ 6 -^-( 1 + (3) (7) 
o + a 

Clearly, if |J| >0, then Tr(J) < 0, and the solution is stable. So we can get 
the critical value A c of A. For simplicity, we let 5 = 1. Then 



(q + l)( 7 + /?) 
AC = (k) {) 

If A < A c , the solution (I) is stable, and the disease will die out. Otherwise 
solution (I) is not stable. 

For solution (II), the Jacobean is: 



(9) 



S\(k)-(5+a)( 1+ p) n 
A(fc)(5+/3) U 

Considering y = ^j^pgy^ > 0, we also can get A c (let 5=1): 

K = M (10) 

When A > A c , the solution (II) is stable, which means that the disease 
will pervade the network; otherwise the disease will disappear. Noticing the 
expressions(8) and (10), we find that A c is a critical parameter. If A < A c , 
the solution (I) is stable, and the disease will disappear from the network; if 
A > A c , the solution (II) is stable, and the disease will spread in the system. 
From (8) and (10), it is obvious that A c is governed by the parameters a, (3, 
7, and (k). We can increase the treatment rate or decrease (k) to raise the 
threshold to prevent disease from spreading. 
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2.2 The spread of disease with treatment on scale-free network 



In the above section, we have studied the epidemic model on regular network. 
But the investigations have shown that a large number of systems, such as in- 
ternet, world- wide- web, physical, biological, and social network, exhibit com- 
plex topological properties [14- 16], particularly scale-free network feature [14]. 
Recent works have examined the spread of computer viruses on the scale free 
networks[7,8,10]. The results show that the intrinsic epidemic threshold is ab- 
sent in both SIS model and SIR model on scale-free(SF) networks. In this 
section, we analyze our modified model on the scale free networks, of which 
the degree distribution is p{k) = C f(k)k~ u , where f(k) is the function of k. 
Suppose x k (t) and y k (t) are the density of the healthy and infected nodes with 
given degree k, and the mean-field equations are[9,10]: 

dx k (t) 

dt = ^(l -x k - y k ) - ax k - \kx k e k (y(t)) + -yy k (11) 
= Xkx k Q k (y(t)) - ( 7 + /3)y k (12) 



where Q k (y(t)) stands for the probability that an edge emanating from a node 
of degree k points to an infected site, Q k (y(t)) = J2 k ' p{k /k)y k '(t), where 
p(k /k) is the probability that a node with k degree points to a node with 
k degree. For uncorrelated networks[17], p{k /k) = kpik )/ (k), which means 
that the probability that a node points to a node with k degree is proportional 
to its degree and the degree distribution p(k ), and (k) is the normalization 
factor. From the definition of Q k (y(t)), we find that it is independent of k for 
uncorrelated networks[17]: 

e k (y(t)) = e(y(t)) = (fc^EfcWWW ( 13 ) 



Let 9X q^ = and 9y ^ = 0, we can get stationary solution: 

Xk = 2+1 (14) 

= AM} 

Vk (i + a )( 7 + /3) + (l + /3)AA;e 1 ' 



In the above expression, we have let 5 — I.. 

Substituting the expression (15) into (13), we get self-consistent equation of 
(-): 

i _ _ \'k 2 e i / \'k 2 e 



= 



(i + P)(k)^ p ( k) i + \'ke (i + p){k) \ i + \'ke) (16) 
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where A' = , u i y, ffl A, 

(l+a)(7+/3) 

We can see that = is a solution of Eq.(16). To allow a nonzero solution 
0(0 G (0, 1]) of Eq.(16), the following inequality must be assumed: 

7 A'A; 2 0) 1 > 1 (17) 



K0- + P) (*> 

From Eq. (17), we get the threshold value of A': 



(k 

where (A; 2 ) = J2kk 2 p(k), then: 



K = (i + P)M- (18) 



A c =(l + ^)(7 + «)^ (19) 



From (19), we can see that A c is dependent on 7, a, (3, and 4§y. If A > A c , 
the disease will spread on the networks, otherwise the disease will die out. We 
now discuss A c for different f(k). 

(I) For f(k) = 5k,k c , then p(k) = Ck~ v 5k,k c (k c > 2). The network is homoge- 
neous, (k) = k C) (k 2 ) = k 2 c) so 

K = ' 1 + ") (7 + /?) (20) 



Clearly, there is a nonzero threshold A c , in agreement with the result on 
the regular network(see Eq. (10)). When A > A c , there is a nonzero = 
X ^\{k\^i+p) +l3 ' 1 °f Eq.(16). A c is an increasing function of 7, a and (3, . We can 
increase the threshold A c by increasing the rate of treatment 7 and decreasing 
k c - 

(II) For f(k) = 1, the network is scale free with a power law degree distribution 
p(k) = Ck- V {v E (2,3]), then 

+00 1 
(k) = k P( k ) - C-m 2 ' v (21) 

k=m ^ 

k 2 ) = k2 P(k) - / k 2 ~ v dk (22) 

k=m Jm 
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As (2 — v) is bigger than — 1, so (A; 2 ) is divergent, and — > 0, then A c — > 
for k — > oo. Therefor the threshold is absent. Thus the treatment is of no 
effect to the disease. 

(Ill) For f(k) = e~ k / kc , f(k) decreases rapidly for k > k c , the network is a 
finite size scale free network [18]. Then 

A c = (l + Q; )(7 + /3) (/// 



V Ll—Vp-k/kc 

(l + a)(7 + /?)-. " ' 



E fc k 2 ~ v e- k / k - 



Where is the incomplete gamma function. The threshold A c is nonzero 

for finite k c . Without surprise, the threshold is an increasing function of 7, a, (3. 
Comparing (23) with (20), one can see, for network with degree distributions 
p(k) = Ck~ v e~ k l kc A c is also very small. From (I)-(III), we find that the 
network degree distribution, in some sense, determines the spread of infectious 
disease. 



3 Conclusion 



To summary, we have suggested an epidemic model with birth rate and death 
rate. Through mean-field analysis, we find that on regular network the epi- 
demic threshold is an increasing function of treatment rate 7, death rates a 
and f3; while for power law degree distribution network epidemic threshold is 
absent in the thermodynamic limit, so that the treatment thus is of no ef- 
fect to the disease, which is the same as the result of standard SIS model. 
So to prevent the infectious disease spreading in the "networks", apart from 
increasing the cure rate, we should pay more attention to the structure of 
"networks". We should point out that we do not analyze specific disease in 
our model. For a specific disease, there is a great need to analyze the disease 
through modelling and comparing the epidemic model with real data. 
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